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The Structure of Particle-Laden Jets 
and Nonevaporating Spraya 


SUMMARY 

The first phase of a study of the structure of sprays, limited to the 
properties of particle-laden jets and nonevaporating sprays, is described. 
Several models of th- se processes were developed. Model performance was 
evaluated by comparison of predictions with existing measurements in 
particle-laden jets, as well as with new measurements in nonevaporating 
sprays completed during this investigation. Analysis and measurements 
were limited to the dilute portions of the flows, where the volume fraction 
of the continuous phase was greater than 99.1%. 

Three models of the process were developed: (1) a locally homogeneous 
flow (LHF) model, where slip between the phases was neglected; (2) a de- 
terministic separated flow (DSF) model, where slip was considered but 
effects of particle/drop dispersion by the turbulence were ignored; and 
(3) a stochastic separated flow (SSF) model, where effects of interphase 
slip and turbulent dispersion were considered using random sampling for 
turbulence properties in conjunction with random-walk computations for 
particle motion. All three models used a k-e model which was extensively 
evaluated for constant and variable density single-phase jets during 
earlier work in this laboratory. 

The new spray experiments employed vacuum pump oil — to insure negligible 
evaporation. The sprays were produced by an air-atomizing injector. Mean 
and fluctuating velocities and Reynolds stress were measured in the con- 
tinuous phase using laser Doppler anemometry. Liquid mass fluxes were 
measured by isokinetic sampling. Drop sizes were measured by slide 
impaction. 

The LHF and DSF models did not provide very satisfactory predictions 
over the present data base. The DSF model generally underestimated the 
rate of spread of the dispersed phase as a result of ignoring effects of 
turbulent dispersion. The LHF model provided reasonably good predictions 
for flows containing tracer-like particles, but was unsatisfactory for 
most practical flows. In earlier evaluations, LHF models generally over- 
estimate the rate of spread of dispersed phases due to neglect of slip. 
Similar behavior was observed during this study for particle-laden jets, 
however, the LHF model underestimated spread rates for the present non- 
evaporating sprays. This indicates that LHF models do not always provide 
an upper bound on the rate of development of dilute particle-laden flows 
as suggested in the past. 

In contrast to the other models, the SSF model provided reasonably 
good predictions over the data base. While this result is encouraging, 
uncertainties in initial conditions for much of the data limits the 
thoroughness of the evaluation of the SSF model. Some effects of 
particles on turbulence properties were observed at high particle mass 
loadings. Treatment of such dense particulate flow effects will require 
extension of the present SSF model. The SSF approach, however, appears 


li 


to provide an attractive formulation for treating nonlinear interphase 
transport processes in particle-li-den turbulent flovs. Current work in 
this laboratory is considering extension of the method to evaporating 
sprays . 
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NOMENCLATURE 


Symbol 

a 
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P 
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G 

P 

k 


L 

e 


m 


P 


m 

o 


M 

o 


n 

n 


i 


Re 
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Description 

acceleration of gravity 

particle concentration 

drag coefficient 

parameters in turbulence model 

injector diameter 

particle diameter 

mixture fraction 

particle mass flux 

turbulence kinetic energy 

dissipation length scale 

particle mass 

injector flow rate 

injector thrust 

number of particles per unit time in group i 

number of particle groups 

Reynolds number 

radial distance 

source term 

particle source term 

time 

eddy lifetime 
particle transit time 
axial velocity 
particle velocity vector 

radial velocity PRECEDING PAGE BLANK NOT FILMED 


Symbol 
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Description 

Favre radial velocity 

tangential velocity 

axial distance 

particle position vector 

path length of particles in an eddy 

time of particle residence in an eddy 

rate of dissipation of turbulence kinetic energy 

turbulent viscosity 

density 

particle relaxation time 
turbulent Prandtl/Schmldt number 


generic property 


Subscripts 

c 

P 

o 

00 


centerline quantity 
particle property 
injector exit condition 
ambient condition 


Superscripts 

( )’ fluctuating quantity 

( ) time mean value 

o 

( )" 


vector quantity 
instantaneous quantity 
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1. INTRODUCTION 


The objective of this investigation was to complete measurements 
of spray structure, useful for evaluation of models of the process. 

The tests considered nonevaporating sprays, generated by an air- 
atomizing Injector, in a still environment. This arrangement has 
simple geometry and well-defined boundary conditions, which facili- 
tates model evaluation. Furthermore, nonevaporating sprays highlight 
effects of drops on the properties of the continuous phase and effects 
of drop dispersion by turbulence in sprays, while minimizing compli- 
cations due to density variations in the flow. The new data was also 
used to begin model evaluation considering methods typical of recent 
spray models. 

Existing data on the mean and turbulent structure of nonevapor- 
ating sprays is limited. A recent review discusses early work in 
the field [1]. Subsequently, Alpert and Mathews [2,3] report pre- 
dictions and measurements for nonevaporating sprays having configur- 
ations encountered in sprinkler systems. Due to the complexity of 
this geometry and the limited measurements available, however, they 
suggest that additional measurements and analysis would be desirable. 
Yule et al. [4] report measurements in nonevaporating sprays from a 
twin-fluid injector in a coflowing stream. Measurements of particle 
size were undertaken using a laser tomographic light-scattering 
technique while mean drop and continuous-phase velocities were mea- 
sured using laror Doppler anemometry (LDA) . These measurements ate 
currently being employed to evaluate spray models in this laboratory. 
The results of this comparison will be reported in the future. 

The present investigation supplements the measurements of Yule 
et al. [4] considering the simpler limiting case of a spray in a 
stagnant environment. Mean and fluctuating velocities of the con- 
tinuous phase were measured, using LDA techniques. Drop sizes were 
measured using the Fraunhofer diffraction and the slide impaction 
methods. Liquid fluxes in the spray were found with an isokinetic 
sampling probe. 

These new measurements, as well as existing measurements of the 
structure of particle -laden jets, were employed to initiate evalu- 
ation of models to predict spray structure. Since these flows only 
Involve Interphase nomen turn exchange, modeling efforts concentrated 
on effects of turbulent fluctuations on momentum exchange between 
the phases as well as the dispersion of the dispersed phase by 
turbulent fluctuations. 

The structure of sprays and other particle-laden flows is 
generally influenced iiy turbulent dispersion of the discrete phase. 
Turbulent dispersion of particles is examined during this investigation 


Numbers in brackets denote references 
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by comparing predictions of several theoretical models with existing 
measurements. The study emphasized solid-particle-laden jets In a 
still environment. These results are also of interest for spray 
modeling, however, since the geometry approximates near injector 
conditions while consideration of solid particles avoids compli- 
cations due to polydlsperse drop size distributions and drop coa- 
lescence. 

Past models of turbulent particle-laden jets often consider 
two limiting cases instead of treating turbulent particle dispersion 
[1]. At one limit the particles and the continuous phase are as- 
sumed to have equal rates of turbulent diffusion. The locally homo- 
geneous flow (LHF) approximation provides a consistent formulation 
of this limit. This implies that interphase transport rates are 
infinitely fast, so that both phases have the same velocity at each 
point In the flow. The LHF approximation provides best results for 
flows containing small particles, where characteristic response 
times cf particles are small in comparison to characteristic times 
of turbulent fluctuations. LHF models have been extensively evalu- 
ated during earlier work in this laboratory,* but only yielded accurate 
predictions for particle sizes smaller than most practical appli- 
cations [5-7]. 

Turbulent particle dispersion is neglected entirely at the 
other limit. This implies that particles follow deterministic 
trajectories since they only Interact with mean properties of the 
continuous phase, yielding a deterministic separated flow (DSF) 
model. Such an approximation is appropriate for flows containing 
large particles, where characteristic particle response times to 
flow disturbances are large in comparison to characteristic turbu- 
lent fluctuation times. Several spray models have been proposed 
along these lines, e.g., El Banhawy and Whitelaw [8], Hongia and 
Smith [9], and Boyson and Swithenbank [10], among others [1]. Due 
to the complexities of sprays and uncertainties in Initial con- 
ditions, however, the value of DSF models has not been clearly 
established as yet. 

Most practical particle-laden flows exhibit properties between 
these limits and require consideration of turbulent particle dis- 
persion. Early dispersion models, discussed by Yuu et al. [11] 
apply a gradient diffusion approximation with empirical correlations 
of turbulent particle exchange coefficients. This approach is not 
practical, however, since such exchange coefficients are influenced 
by both particle and turbulence properties — requiring excessive 
effort to accumulate a data base sufficient for general appli- 
cation of the method. 

Several recent studies of turbulent particle dispersion use 
stochastic separated flow (SSF) methods to circumvent the limi- 
tations of the gradient diffusion approach. Stochastic analysis 
requires an estimate of the mean and turbulent properties of the 
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continuous phase. Particle trajectories are then computed by random 
sampling to find Instantaneous continuous phase properties. Mean 
and fluctuating particle properties are found by Monte Carlo aethoda — 
where a statistically significant number of particle trajectories 
are averaged to obtein system properties. 

SSF models have been applied to particle-laden jets. Yuu et al. 
[11] use empirical correlations of mean and turbulent properties 
for SSF analysis of their pertide dispersion measurements. Gosman 
and Ioannldes [12] propose a more comprehensive approach, where flow 
properties for the stochastic calculations sre computed with a k-e 
turbulence model. This approsch is examined here, after only minor 
modification. 

In the following, the models are described first of all, fol- 
lowed by initial evaluation using existing measurements for particle- 
laden jets. Experimental methods used for measurements in nonevapor- 
ating sprays are then discussed. The report concludes with an evalu- 
ation of the models using the spray measurements. Other reports of 
the findings of this investigation can be found in Refs. [13-15]. 


2. THEORY 


2.1 General Description 

A k-e turbulence model was used to find the properties of 
the continuous phase for all models, since this approach has been 
thorough 'y calibrated for both constant and variable density jets 
[5-7,16]. All conditions to be considered can be modeled as steady 
axlsymmetric boundary layer flows. Nozzle exit Mach numbers were 
less than 0.3; therefore, kinetic energy and viscous dissipation of 
the mean flow, as well as gas density variations, were neglected with 
little error. Other assumptions will be described when introduced 
since they differ for each model. 

2.2 LHF Model 


The formulation of the LHF model corresponds to the general 
treatment of the continuous phase for all three models. The LHF 
approximation Implies that both phases have the same velocity at 
each point in the flow, i.e., local thermodynamic equilibrium is 
maintained. Therefore, the flow corresponds to a variable density 
single-phase fluid due to changes in particle concentration even 
though the density of each phase is constant. Following past work 
[5-7,16] the governing equations are written in Reynolds averaged 
form, as follows: 

3 /"ax ^13/ — o . . 1 9 , P t 3$. a * 

a? (pu<p) + 7 3? (rpv ♦> “ 7 37 (r 57 a? + s * (1) 

9 


4 


where 
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Pv° ■ pv + p'v' (2) 


The peremeters $ and appearing in Eqs. (1) and (2) are summarised 
in Table 1, along with the appropriate empirical constants. The 
empirical constants were established for single-phase flows [5-7,16]. 
They were not changed during present calculations for particle-laden 
jets. The turbulent viscosity was calculated from the turbulent 
kinetic energy and its rate of dissipation as follows: 

M t " C y Pk 2 /e (3) 


For present assumptions, the instantaneous particle concentration 
and flow density are only functions of mixture fraction — corresponding 
to the equilibrium state reached when f and (1-f) kg of Injected and 
ambient fluid, respectively, are adiabatically mixed. This yields 


C/C q - f ; p " 1 - f/p o + (l-f)/ Poo (4) 

Since C and p * are linear functions of f in the domain 0 £ f 1, 
their mean values can be found by substituting f in Eq. (4). In this 
instance, it is not necessary to solve a transport equation for mixture 
fraction fluctuations and to adopt a probability density function (PDF) 
for f — which is usually necessary during LHF analysis of sprays [5-7]. 

Ambient values of u, f, k and e are all zero for the flows treated 
here. Gradients of these quantities are also zero at the axis. 

Initial conditions varied over the data base. For the round jets 
and the spray, initial conditions were approximated as either slug 
flow or fully developed pipe flow. 

When slug flow is assumed at the injector exit, initial conditions 
were prescribed similar to past work [5-7]. All properties were 
assumed to be constant, except for a shear layer having a thickness 
equal to 1 percent of the injector radius at the passage wall. The 
constant property portion of the flow was specified as follows: 


x - 0 , r < 0.99d/2; 

u“M/m , f ■ 1 , k ■ (0.02u ) ^ , 

o o o o o o 

e = 2.84 x 1C)" 5 u 3 /d 

o o 


( 5 ) 
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Table 1 Source Terms in Eq. (1) 


$ 


S 


♦ 


+ * (p - P) 


<£> 2 - * 


% »« I # 2 


C e 2 P k 


Notes: 

1. Positive sign is used In S- for vertical upward 

flow. u 

2. Turbulence model constants are assigned the 
following values: 

- 0.09, C £ - 1.44, o fc - 1.0, 

a t ■ 1.3, Oj"0.7, and 

C * 1 . £ ! 4 (variable density flows) ■ 1.89 
e 2 (constant density flows). 
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Equation (5) provides the inner boundary conditions of the shear layer 
until it reaches the jet axis. The initial variation of u and F is 
taken to be linear in the shear layer. Initial values of k and e in 
the shear layer were found by solving their transport equations while 
neglecting convection and diffusion terns. 

When fully-developed pipe flow was assumed at the injector exit, 
f was taken to be unity while u was obtained fron the power law 
expressions provided by Schlictl&g {17] — allowing for variation of 
the power with Reynolds number. Initial values of k and e were 
obtained from Hinze [18] for fully-developed pipe flow in tRe present 
Reynolds number range. 

The calculations were performed using a modified version of 
GENMIX [19]. The computational grid was similar to past work [5-7]: 

33 cross-stream grid nodes; streamwlse step sizes limited to 6X of 
the current flow width or an entrainment increase of 10Z— whichever 
was smaller. 

2.3 DSF Mo del 

Continuous Phase . Both separated flow models adopt the 
main features^of the LHF model, but vnly for the gas phase. The 
particle/drop phase is treated by solving Lagrangian equations of 
motion for the particles and then computing the source terms which 
appear in the governing equations for the gas phase due to inter- 
phase transport processes. This procedure corresponds to the 
particle-tracking or particle-source-in-cell methods used in most 
recent models of sprays [1-3,8-10]. With this approach, particles 
leaving the jet exit are divided into n groups, defined by initial 
position, velocity, and direction. 

The void fraction of the flows was always greater than 99. IX; 
therefore, the volume occupied by the discrete phase was ignored 
with little error. Under these conditions it it also reasonable to 
neglect the effect of particle motion on turbulence properties — 
although exchange of mean momentum between phases was considered as 
noted earlier. The gas-ghase density is constant, simplifying Eq. 

(1), and a solution for f is no longer needed since particle con- 
centrations are found from the discrete phase solution. 

The interaction between the dlscrete_and continuous phases yields 
an additional source term in Eq. (1) for u. This term is found by 
computing the net change in momentum of each particle group 1 as it 
passes through a computational cell j (only one index is needed to 
define a cell since the calculation is parabolic and each cell is 
defined by its radial position). This yields the following source 
term expression 


* 

In the following, parvicle will denote drops (assumed to be 
spherical) as well as solid particles. 
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s n • k ■ 


pi (^ U pi^in “ ^ U pi^out)j 


where n . is the number of particles per unit time in each group, and 
"in" ana "out” denote conditions entering and leaving a computational 
cell. 


The continuous-phase equations were solved in the same manner as 
the LHF model. 


Discrete Phase . The main assumptions of the particle trajectory 
calculations are as . ollows: due to the high void fractions of the 
measurements, dilute-particle flow was assumed with drag equivalent 
to a single particle in an unbounded environment and negligible 
particle collision or coalescence when drops are considered; drag was 
treated empirically, assuming quasisteady flow for spherical particles; 
since p / p > 200 for conditions treated here, effects of virtual mass 
Bassett p forces, Magnus forces, etc., can be neglected with little 
error; effects of turbulent fluctuations were neglected, similar to 
most separated flow models; and ambient conditions of the particles 
were taken to be local mean flow properties. The last two assumptions 
are characteristic of the DSF formulation 11-3,8-10] and will be 
relaxed for the SSF model. Tha remaining assumptions a:, a typical of 
separated flow models of dilute sprays and are disucsted more com- 
pletely elsewhere [1]. 

After adopting these assumptions, the pcsiti^n and velocity of 
each particle group is found by integrating 


dx 

—El . 

dt 


Pi 


i - 1,2,3 


(7) 


dt 


3 PC D 

4 d P 

P P 


(,J i " u p i ) l“ ~ “ p l + *i • i ” 1 * 2 * 3 


( 8 ) 


Particle Reynolds numbers did not reach the supercritical flow 
regime; therefore, the standard drag coefficient for solid spheres 
was approximated as follows [1]: 


s, 


) 


Re < 1000 


( 9 ) 


0.44 


Re > 1000 
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The particle motion equations, Eqs. (7) and (8), mere solved at 
the same time as the gas-phase equations, in a stepwise fashion. A 
second-order finite difference algorithm was employed for these com- 
putations. The computations employed no less than 100 (particle- 
laden jets) and 400 (sprays) groups. 

2.4 SSF Model 


The SSF model involves finding trajectories of a statis- 
tically significant sample of individual particles as they move away 
from the injector and encounter a random distribution of turbulent 
eddies — using Monte Carlo methods. The treatment of the continuous 
phase was identical to the DSF model. 

Key elements in the SSF model are the method used to specify 
eddy properties and the time of interaction of a particle and a 
particular eddy. The approach used to find these properties follows 
Gosman and Ioannides [12], but differs in some details. Properties 
were assumed to be uniform within each eddy and to change randomly 
from one eddy to the next. Particle trajectory computations were the 
same as the DSF model, involving solution of Eqs. (7) and (8); how- 
ever, mean gas properties in these equations were replaced by the 
instantaneous properties of each eddy. 

The properties of each eddy were found at the start of particle/ 
eddy Interaction by making a random selection from the probability 
density function (PDF) of velocity. Velocity fluctuations were 
assumed to be isotropic with a Gaussian PDF having a standard devi- 
ation of (2k/3)^'^ and mean components u, v, 0. The cumulative 
distribution function for each velocity component was constructed 
and sampled. This involved randomly selecting three numbers in the 
range 0-1 and computing tne velocity components for these three values 
of the cumulative distribution function. 

A particle was assumed to interact with an eddy for a time which 
is the minimum of either the eddy lifetime or the transit time re- 
quired for the particle to cross the eddy. These times were esti- 
mated assuming that the characteristic size of an eddy is the dissi- 
pation length scale. 


L 

e 



k 3 ' 2 /e 


( 10 ) 


Gosman and Ioannides [12] compute the eddy lifetime as t * L /|u’|, 
however, v>e obtained better "greement with measurements using 


t - L / (2k/3) 1/2 
e e 


(id 
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A portion of th« piaitnt calculations employed tha method propoaad by 
Goamen and Ioannidas to find tha transit time. In this csss, tha 
transit tins of a partlda is found from tha linsarissd aquation of 
partlda notion in a unifora flow (LSSF nodal) 

t t - -t In (1 - L e /(ifu" - u"|)) (12) 

-a 

where u" - u" is tha velocity diffaranca at tha start of tha Intar- 
act ion and p is the llnaarited particle relaxation tins 

T " 4 pd p /<3 C D |> - Sjl) (13) 

When L > t|u" - u"| , the linearised stopping distance of tha particle 
is saafler than thl characteristic length scale of tha eddy and Eq. 
(13) has no solution. In this case, tha eddy has captured tha par- 
ticle and the interaction tine is the eddy lifetime. 

The linearization of the particle motion aquation is an unde- 
sirable feature of LSSF model, since nonlinear effects are considered 
elsewhere in the analysis. Therefore, a second method waj also used 
where particles were assumed to Interact with an eddy as logg as the 
time of interaction, At, and the distance of interaction |Ax J, 
satisfied the following criteria (NSSF model) p 

At < t , |Ax I < L (14) 

— e 1 p 1 — e 

Particle capture by an eddy corresponds to ending the interaction 
with the first criterion while a particle traverses an eddy when the 
interaction is ended with the second criterion. 

The remaining computations are similar to the DSF model, except 
that more particle trajectories must be considered to obtain statis- 
tically significant particle properties (generally 1000-2500 tra- 
jectories were used). A by-product of the additional calculations, 
however, is that the SSF models yield both mean and fluctuating 
particle properties. This provides an additional test of model 
performance. 


3 . MODEL EVALU ATIO N: PARTICLE-LADEN JETS 
3 . 1 D ata Ba se 

The data base used to develop and evaluate the particle-laden 
Jet models is summarized in Table 2. The findings of Hinse (18] and 
Snyder and Lumley [20] were used for Initial evaluation of the SSF 
model — similar to Gosman and Ioannides [12]. The other measurements, 
for particle-laden Jets, were used for model evaluation. 


Table 2 Summary of Measurements Used for Evaluation 
in Solid Particle-Laden Flows 3 
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Continuous phase was air at normal temperature and pressure. Round jets were injected into still air. 
Ratio of the mass flow rates of particles and air. 

Slip only considered for separated flow predictions. Slip conditions were corrected to an equivalent 
initial momentum and no slip for LHF model predictions. 
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Adequate Information concerning initial conditions was not 
available for any of the particle-laden jet experiments. Present 
estimates of initial conditions for the computations are summarized 
in Table 2. The jet was formed by relatively short nozzles for the 
measurements of Yuu et al. [11] and Me Comb and Salih [21,22]; there- 
fore, the slug flow approximation was employed to specify these 
initial velocity profiles. Laats and Frishman [23,24] and Levy and 
Lockwood [25] used relatively long constant diameter tubes for in- 
jectors and fully developed pipe flow was taken to be the best ap- 
proximation for properties at the jet exit. Due to lack of other 
information, particle concentrations were assumed to be uniform 
over the jet exit in most cases. The data of Laats and Frishman 
[23,24] was an exception, they provided the ratio of centerline to 
average particle mass flux at the jet exit; therefore, simple power 
law distributions were assumed which yielded this ratio. Yuu et al. 
[11] and Levy and Lockwood [25] provide sufficient information con- 
cerning properties of their particle mixing and injection system so 
that Initial particle velocities could be estimated. A DSF model 
was used for these calculations. For lack of such details, it was 
necessary to assume no slip between the phases at the injector exit 
for the measurements of Refs. [21-24]. Effects of initial particle 
velocity fluctuations were ignored in the separated flow calculations, 
but are intrinsically considered in the specification of initial con- 
ditions for LHF calculations. 

3.2 Model Calibrations 


The evaluation of the k-e turbulence model predictions for 
constant and variable density single-phase jets is described else- 
where [1,5-7,16]. The data base used for this evaluation included 
measurements by Shearer et al. [5], Mao et al. [6], Wygnanskl and 
Fiedler [26], Becker et al. [27], Corrsin and Uberoi [28], and 
Hetsroni and Sokolov [29] . The model generally yielded satisfactory 
predictions of u, F, k and Reynolds stresses in these flows. There- 
fore, its use to find mean and turbulent properties of the continuous 
phase during present calculations is justified, although further im- 
provements could be adopted by the LHF, DSF and SSF models with no 
fundamental difficulty. 

The method for specifying eddy properties and lifetime in the 
SSF model is somewhat arbitrary. Therefore, the procedure was evalu- 
ated using the analytical results of Hinze [18] for dispersion of 
infinitely small particles injected at a constant rate from a point 
source in a homogeneous isotropic turbulent flow — similar to Gosman 
and loannides [12]. The calibration results are illustrated in Fig. 1. 
Stochastic model predictions from Gosman and loannides [12] and the 
present model are shown, along with the exact analytical results from 
Hinze [18]. The present model is in good agreement with the exact 
results, indicating satisfactory calibration. The results for the 
LSSF and NSSF models are the same in this case, since infinitely small 
particles are always captured by eddies and only the method for 


C/C 


12 


0RKM4M- * 

rm POOR QURUIt 


1 1 1 1 1 1 I 



r(mm) 


Fig. 1. Analytical and stochastic solutions for 
the dispersion of small particles in 
homogeneous isotropic turbulent flow 
after long diffusion times. 
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computing eddy lifetime is relevant. The Gosman and Ioannides [12] model, 
as posed, yields poorer agreement due to a computational error in 
their original work that has since been corrected [30]. 

Predictions of the present stochastic LSSF model, as well as the 
version of Gosman and Ioannides [12] and the measurements of Snyder 
and Lumley [20] °re illustrated in Fig. 2. These experiments involved 
dispersion of individual parti - les which were isokinetically injected 
into a uniform turbulent flow downstream of a grid. The two models 
yield nearly the same results for heavy particles where the cali- 
bration error of [12] has little effect. The present model yields 
better agreement for the lighter particles where calibration of 
characteristic eddy time and length scales is more significant since 
light particles tend to be captured by the eddies. Results for the 
NSSF model were also compared with this data, these predictions were 
slightly improved over the LSSF results illustrated in Fig. 2. 

3.3 Results and Discussion 


Measurements of Yuu et al. [11]. The significance of tur- 
bulent particle dispersion in particle-laden jets can be appreciated 
from the results appearing in Fig. 3. The particle concentration mea- 
surements of Yuu et al. [11] are illustrated along with the predictions 
of the LHF, DSF and NSSF models (NSSF and LSSF model predictions are 
nearly identical in this case) . Predictions are provided for values 
of x/d which bound the range of the data. 

The rate of particle spread is overestimated by the LHF model 
in Fig. 3, since the particles are predicted to be too responsive to 
lateral turbulent fluctuations as a result of neglecting slip between 
the phases. On the other hand, the rate of particle spread is under- 
estimated by the DSF model, since radial particle velocities which 
contribute to particle spread are only generated by the v velocity 
component — which is small in comparison to radial velocity fluctu- 
ations near the axis of jets [26]. The DSF model predictions also 
indicate that particles are confined to a progressively narrower 
range of r/x as x/d increases, which is contrary to the measurements. 
The DSF model yielded poor results for the remainder of the evaluation 
and will not be considered in subsequent figures for particle-laden 
jets. In contrast to the relatively poor performance of the LHF and 
DSF models, the SSF model predictions are in good agreement with the 
measurements in Fig. 3. 

Predictions of the LHF and LSSF models are compared with data 
from Yuu et al. [11], over a broader range of conditions, in Fig. 4 
(the results for u. 50 m/s can be compared with Fig. 3 to indicate 
the differences between LSSF and NSSF predictions). In general, the 
LSSF model provides good predictions while the LHF model continues to 
overestimate the rate of particle spread. 

Measurements of McComb and Salih [21,22]. The LHF and SSF (both 
SSF versions yielding nearly the same results) model predictions are 



Predicted and measured particle dispersion in a uniform grid-generated 
turbulent flow. Data of Snyder and Lumley [20], 





15 


original page is 

OF POOR QUALITY 



r/t 


Fig. 3. Comparison of LHF, DSF and NSSF model 

predictions of particle dispersion in an 
axisymmetric jet — data of Yuu et al. [11]. 
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compared with the measurements of McComb and Salih [21 9 22] in Figs. 

5 and 6. Results for relatively small titanium dioxide particles, 

2.3 Mm in diameter , appear in Fig. 5. In this case, these small 
particles act nearly like tracer particles and the LHF model yields 
adequate results. The SSF models underestimate particle dispersion 
slightly, but are also in reasonably good agreement with the mea- 
surements. 

Results for larger particles, 5.7 pa in diameter, are illus- 
trated in Fig. 6. These particles were tungsten and had a density 
roughly five times greater than the titanium dioxide particles con- 
sidered in Fig. 5, cf. Table 2. The greater size and density of the 
tungsten particles increases their inertia so that effects of slip 
are more important for the conditions of Fig. 6. In this case, only 
the SSF models yield adequate agreement with the measurements while 
the LHF model overestimates both the rate of spread of the particles 
and the flow. 

Results thus far were limited to relatively low particle load- 
ings. This implies that while the gas flow influences particle 
dispersion, the effect of the particles on the structure of the 
gas flow was small. This deficiency in the evaluation of the 
models is rectified by comparison with the measurements of Lasts 
and Frishman [23,24] and Levy and Lockwood [25] in the following. 

These measurements involved relatively large particle loadings, re- 
sulting in significant effects of particles on the structure of the 
gas phase. 

Measurements of Lasts and Frishman [23,24]* Particle mass 
fluxes along the axis were observed to Increase for a time near the 
jet exit for a portion of the data of Lasts and Frishma r [23,24]. 

The authors attribute this effect to Magnus forces Induced by the 
particle/gas mixing and injection processes. This effect could not 
be prescribed sufficiently to include it in the calculations; there- 
fore, data exhibiting such trends are not considered in the following. 

Predictions of the NSSF and LSSF models are compared with the 
measurements in Figs. 7 and 8. Only a range of initial gas velocities 
were specified for the data; therefore, results are provided for both 
limits of the velocity range. The predicted effects of initial ve- 
locity changes in this range are small, in agreement with observa- 
tions [23,24]. 

Predicted and measured mean gas velocities along the jet axis are 
illustrated in Fig. 7. Results are shown for both particle-laden jets 
and an air jet. The gas is not strongly Influenced by the momentum of 
the particles at low particle loadings and mean-gas velocities approach 
the air-jet properties — even for large particles. As particle loading 
increases, however, the rate of decay of centerline velocity is reduced. 
This trend is reproduced by the NSSF model, with the LSSF version 
yielding essentially the same results — except at the highest loading 
ratios where the models overestimate the rate of flow development. 
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Fig. 6. Predicted and measured axial gas velocities and 
particle concentrations in particle-laden Jets — 
data of McComb and Salih [21,22]. 
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Predicted and measured particle mesa flux for various conditions 
are illustrated in Fig. 8. The results are for x/d • 28. 5 since Lasts 
and Frishmsn suggest that this is beyond the range where Magnus effects 
were important. The LHF model significantly overestimates the rate of 
particle dispersion in all cases since the test particles have signifi- 
cant inertia. The SSF models, however, yield reasonably good pre- 
dictions — particularly the NSSF version — except at the highest particle 
loadings where the rats of development of the flow is overestimated 
similar to Fig. 7. Less satisfactory agreement is also observed for 
the 80 pm diameter particles at a loading ratio of 0.3, where pre- 
dictions underestimate the rate of particle spread. In this case. 
Initial slip may be a factor since such large particles are less likely 
to equilibrate in the injector tube, e.g., use of 30% slip at the in- 
jector exit yields good agreement between predictions and measurements. 

Two effects may be responsible for the poorer performance of the 
SSF models at higher loading ratios, seen in ELgs. 7 and 8. The effect 
of Magnus forces, discussed earlier, causes particles to be deflected 
toward the centerline of the flow in a manner not considered in the 
analysis. This would reduce rates of apparent turbulent particle dis- 
persion and would also feed additional momentum from the deflected 
particles to the flow near the centerline — tending to reduce rates of 
axial velocity decay. A second effect could be the turbulence modu- 
lation phenomena discussed by A1 Taweel and La'tdau [31]. This effect 
involves damping of turbulent fluctuations by the drag of the particles, 
tower turbulence levels reduce rates of particle spread by turbulent 
dispersion as well as the turbulent mixing of the continuous phase, 
which would also explain why the SSF models overestimate rates of flow 
development at high loading ratios. Unfortunately, Information avail- 
able is not sufficient to evaluate the relative significance of these 
two effects. The unusual properties of the measurements of Refs. 

[23,24] even at moderate particle loadings, however, cannot be attri- 
buted to turbulence modulation. This suggests that Magnus forces may 
be a major contributor to the difficulty. 

Measurements of Levy and Lockwood [25]. Predictions of the NSSF 
model (LSSF model predictions are nearly the same) and the measure- 
ments of Levy and Lockwood [25] are illustrated in Figs. 9 and 10. 

LHF model predictions for these flows overestimated rates of flow de- 
velopment similar to results discussed earlier and are not shown on 
the figures. Predictions of fluctuating gas velocities were obtained 
using the assumption ot isotropic turbulence u'^ ■ 2k/3. 

Both predi. tions and measurements indicate relatively small ef- 
fects of particle size and loading on gas flow properties in Fig. 9. 

This is due to the fact that large particle loadings correspond to 
large particles for these data. Therefore, momentum exchange between 
the phases was relatively small for the limited axial distance con- 
sidered during these experiments. Predictions are in reasonably good 
agreement with measurements, except for the highest particle loading, 
where measured turbulence intensities tre underestimated by the 
theory. 
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Fig. 9. Predicted and measured mean and fluctuating 
gas velocities in particle-laden jets — 
data of Levy and Lockwood [ 25 ] . 
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Predicted and measured mean and fluctuating particle velocities 
are illustrated in Fig. 10. Predictions of mean particle velocities 
are in good agreement with the measurements at all conditions. Pre- 
dictions of particle fluctuations, however, are underestimated except 
for the smallest particles in the data base. This discrepancy in- 
creases with the size of the particles, but is not influenced very 
much by the particle loading ratio. 

The progressive underestimation of particle fluctuations with 
increasing particle size can best be explained as artifact of the 
particle/gas mixing and injection system (which was a worm-gear drive 
followed by a relatively short length of ducting). During these 
tests, the large particles have a relatively small velocity change 
as they pass from the jet exit to the measuring station, so that they 
would tend to maintain fluctuation levels which are similar to those 
at the jet exit. In contrast, the small particles interact with the 
flow field to a greater degree, so that effects of the injection 
system are damped out. Since present calculations ignore initial 
particle fluctuations, results are most satisfactory for small par- 
ticles where effects of initial conditions are less persistent. Par- 
ticle velocity fluctuations at the jet exit could also be responsible 
for the somewhat larger velocity fluctuation levels observed for 
large particles in Fig. 9. 

Another explanation for the larger particle and gas velocity 
fluctuations observed for large particles could be that the larger 
slip of these particles enhances turbulent fluctuations in a manner 
that is not considered in the model. Additional data, including 
complete specification of particle and gas properties at the jet exit, 
are needed to clarify the relative importance of turbulence enhance- 
ment/modulation and initial condition effects for these data. 

In fact, the present evaluation of the DSF and SSF models was 
inhibited throughout by insufficient information concerning initial 
conditions. Until these initial condition difficulties are resolved 
by new experiments, decisive evaluation of turbulent particle dis- 
persion models cannot be completed. 


4 . SPRAY MEASUREMENTS 


4 . 1 Test Apparatus 

Figure 11 is a sketch of the apparatus and major instru- 
mentation for the nonevaporating spray experiments. An air-atomizing 
injector, cf. Table 3 for specifications, was mounted at the top of 
a screened enclosure (1 m square and 2.3m high) to reduce room dis- 
turbances. All optical instrumentation was mounted rigidly. There- 
fore, the injector was traversed to obtain profiles of flow quantities 
(two directions using the injector mount a^d the third major traverse 
by moving the entire screened enclosure). The inlet of the exhaust 
system was screened and located 1 m below the plane of the measure- 
ments. Tests showed that operation of the exhaust system had negli- 
gible effect on flow properties at the measuring locations. 
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Table 3 Summary of Test Conditions 
for Nonevaporating Sprays 


Injected fluid 
Injector flow rates, mg/s 


Air Jet 


air 


Case 1 


Sprays 


Case 2 


air and oil** air and oil^ 


Gas 

338 

338 

216 

Liquid 

0 

600 

1400 

G 

Loading ratio 

0 

1.78 

6.48 

Jet momentum, mN 

120.1 

137.2 

70.1 

Initial velocity, m/s^ 

175 

146 

43.4 

Reynolds number** 

2.6 x 10 4 

3.0 x 10 4 

2.4 x 10‘ 

SMD, )jm e 

0 

30 

87 

Spray angle^ 

_ 

30 

33 


All flows employ Spraying Systems alr-atomlzing injector (Model 
1/4J2050 nozzle. No. 67147 air nozzle, 1.194 mm Injector exit 
diameter) . Ambient and injector inlet temperature 300 + 1 K, 
ambient pressure, 97 kPa. 

^Sargent -Welch Scientic Co., Duo Seal Oil, Cat. No. 1407K25, 
density * 878 kg/m 3 , vapor pressure at 38°C » 4 x 10~* mm Hg. 

£ 

Ratio of injected liquid to gas flow rates. 

d Assuming LHF . The viscosity of air was employed for the Reynolds 
number. 

"Measured with the Malvern, Model 3300 Particle Sizer at x/d * 12.6. 
^Determined from liquid flux measurements, cf. Section 4.2. 
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The injector flow system is illustrsted in Fig. 12. Filtered 
dry air is supplied from a storage tank. The flow rate of air is 
controlled with a pressure regulator and metered with a critical- 
flow orifice to insure stable long-term operation. The liquid is 
stored in a tank under pressure, however, the tank is not agitated 
and pressure levels are moderate (0.3-0. 8 MPa). Therefore, the 
dissolved air content of the liquid is negligible. The flow rate 
of liquid is controlled with a valve and metered with a rotameter. 
In order to maintain repeatable flow and atomization conditions the 
entire test cell is heated above normal ambient temperatures — to 
27 +1°C. 

The spray liquid was vacuum-pump oil. This precluded signifi- 
cant evaporation in the flow field — even when finely atomized — due 
to low liquid volatility, cf. Table 3. 

4.2 Instrumentation 


A single-channel, frequency-shifted LDA was used for mea- 
surements of mean and fluctuating gas velocities. Several beam 
orientations provided measurements of all velocity components as 
well as the Reynolds stress. Concentration biasing and effects of 
drops were avoided by employing high concentrations of seeding 
particles (0.2 ym diameter), similar to past work [5]. Data 
densities were high, allowing analog processing to obtain time- 
averaged quantities. 

Injector operation was monitored continuously by measuring drop 
sizes with a Malvern Model 3300 particle sizer. This instrument 
operates on the principle of Fraunhofer diffraction of laser light. 

The measuring region of the instrument included the entire spray 
width, centered at x/d * 12.6, with the input laser beam having a 
diameter of 9 mm. Injector properties were also monitored by 
periodically measuring injector thrust with an impact plate [5]. Spray 
angles were found from liquid flux measurements (taken as the angle 
from the injector exit which bounds the region where the liquid mass 
flux is greater than 1% of its centerline value at x/d * 50). 

Profiles of drop-size distributions were measured by a slide 
impactor, developed earlier [5]. The samples measured on the slide 
were corrected for the collection efficiency of the probe, using a 
correlation due to Ranz, cf. Ref. [5] and references cited therein 
for a description of this correction. Depending on the spray 
condition, 500-3000 drops were measured to establish statistically 
significant drop-size distributions. 

Liquid fluxes were measured with an isokinetic sampling probe. 

The probe had a 3 mm ID sampling port, followed by a gas-flushed 
diverging section which prevented premature drop impaction on probe 
surfaces. The drops were captured on a filter, which was weighed 
after a timed period of collection to determine liquid fluxes. 

Mass conservation checks of liquid flow rates were satisfied within 
15% for the measurements reported here. 
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4.3 Test Conditions 


Three test conditions were considered: (1) e pure air 

jet, formed by the injector, to serve as a baseline; (2) a spray 
having a nominal Sauter mean diameter (SMD) of 30 ym, similar to 
the spray considered in Ref. [5]; and (3) a spray having a SMD of 
87 ym, which exhibits significant effects of slip. The properties 
of these three test conditions are summarized in Table 3. 


5. MODEL EVALUATION: SPRAYS 


5.1 Air Jet Calibration 


The measurements of mean and fluctuating velocities in the 
air jet were in good agreement with earlier measurements by Shearer 
et al. [5] using a similar twin-fluid Injector. The comparison be- 
tween predictions and measurements was also satisfactory — similar 
to Ref. [5]. This established an acceptable baseline for measurements 
in the sprays and these results will not be considered any further here. 

5.2 Initial Conditions 


Due to its small size, measurements at the exit of the in- 
jector were limited to mass flow rate and momentum of the two-phase 
jet. The latter measurement was completed using an Impact plate, 
similar to Ref. [5]. 

Photographs of the two sprays appear in Figs. 13 and 14. De- 
termination of spray angles from photographs is a subjective matter, 
also influenced by the lighting, film and exposure of the film. The 
values fougd from present photographs for case 1 and 2 sprays are 
24 and 21 . These values are less than the angles found from the 
liquid flux measurements — which provide far greater sensitivity for 
determining the presence of liquid. 

Due to limitations of spatial resolution, profiles of spray 
properties were only measured at x/d ^50. In particular, detailed 
measurements of drop-size distributions , liquid flux, and mean and 
fluctuating gas velocities were undertaken at x/d » 50, in order to 
define initial conditions for the SF computations. 

Drop size distributions for various radial positions at x/d ■ 50 
for the two sprays are illustrated in Fig. 15. Although the distri- 
butions differ in detail at the various radial locations, the SMD is 
relatively uniform across both sprays, cf. Fig. 15. The values of 
the SMD at this axial station are larger than the nominal size ob- 
tained by the Fraunhofer diffraction measurements, since: (1) this 
position is farther from the injector which allows drop sizes to in- 
crease by collisions and coalescence in the dense portion of the 
spray; (2) the slide impactor has reduced collection efficiencies for 
small drops [5], which tends to bias its measurement toward larger 
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Fig. 13. Photograph of the case 1 spray. 
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Fig. 15. 


Drop size distributions for the 
nonevaporating sprays at x/d ■ 50 
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drop sices even after correction; and (3) tha slide lmpactor la more 
effective for finding occasional large drops than tha Praunhofar 
diffraction method, which also tends to Increase tha SMD measured by 
impaction. 

The SMD measurements illustrated In Fig. 16 Indicate that drops 
are observed for r/x > 0.2, which Is beyond tha edge cf «>oat fully- 
developed, single-phase jets. This effect is attributable to both 
Initial radial velocities produced at the Injector and dispersion of 
drops by turbulent fluctuations. Further Implications of this ob- 
servation and other measurements at the Initial condition (x/d ■ 50) 
will be considered later. 

5 . 3 Results and Discussion 


Once initial conditions were established, the remaining 
measurements were used to evaluate model predictions. Since the LHF 
model was not limited to a dilute spray and had appropriate initial 
conditions at the injector exit, results for this model were obtained 
throughout the entire flow. The DSF and SSF model predictions, how- 
ever, were limited to the dilute portion of the spray — beginning at 
x/d ■ 50 where initial conditions for these models were measured. 

Predicted and measured mean gas velocities along the axis of the 
two sprays are Illustrated in Fig. 17. The LHF predictions of mean 
velocity for the more finely atomized spray tend to underestimate the 
measurements slightly, similar to the findings of Shearer et al. [5] 
for f.n evaporating spray having a similar SMD. It will become evident, 
however, that this agreement is partly fortuitous, since the LHF 
model significantly underestimates the width of this spray. Typical 
of earlier experience with LHF models of sprays [5] the LHF model over- 
estimates the decay of axial velocity in the coarser spray, due to the 
neglect of effects of slip between the phases. In contrast, the SSF 
model provides satisfactory predictions of mean centerline velocities 
for both sprays. 

The same models are compared with measurements of liquid mass 
flux along the axis of the two sprays In Fig. 18. The LHF models do 
not approach the measurements until large values of x/d, where effects 
of interphase slip are relatively small. The SSF model, however, 
provides reasonably good predict tons throughout the range of x/d 
where it was used. 

Predicted and measured radial profiles of mean axial gas veloci- 
ties are illustrated in Fig. 19, including the initial condition, 
x/d ■ 50. Predictions from both the LHF and SSF models are shown on 
the plot. The results illustrate an interesting property of these 
sprays. During the computations for particle-laden jets, the LHF 
model invariably overestimated the width of the flow since neglecting 
slip causes the rate of dispersion of heavy particles by the turbu- 
lence to be overestimated. In the present case, however, the spray 
spreads more rapidly than the LHF prediction, due to enhanced dis- 
persion of drops. The increased response of the drops in the present 
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Fig. 18. Predicted and measured mean liquid flux along 
the axis of the nonevaporating sprays. 
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sprays , in comparison to the solid particles in the particle-laden 
jets can be attributed to: (1) the smaller density of the drops; (2) 
the greater rates of flow deceleration due to the smaller injector 
dimensions of the sprays; and (3) the larger initial slip between 
the liquid and gas phases in the sprays. It is very encouraging 
that the SSF model can reproduce this unexpected effect, with no 
change in the modeling procedure and empirical constants. 

Predicted and measured radial profiles of liquid m«ss flux are 
illustrated in Fig. 20. In this case, results for the DSF model are 
shovr* along with the LHF and SSF predictions. The DSF model yields 
poor results, similar to particle-laden jets. Neglecting drop dis- 
persion by the turbulence causes the rate of spread of the flow to 
1 > substantially underestimated — even after allowing for the apparent 
radial velocity of the drops at x/d = 50. Other predictions of the 
DSF model were also not very satisfactory; therefore, this method 
will not be considered any further here. 

As before, the SSF model yields the most satisfactory predictions 
of flow properties for the results illustrated in Fig. 20. The per- 
formance of the SSF model, however, is poorer for liquid flux than for 
other measurements considered during this study. Since liquid flux 
predictions are more sensitive to estimations of initial conditions 
than other spray properties, uncertainties in initial conditions are 
a potential source for these errors. Additional measurements of drop 
sizes and velocity will be required to resolve this effect. 

All three components of the velocity fluctuations were measured, 
allowing k to be computed for comparison with predictions. These re- 
sults are illustrated in Fig. 21. The agreement between predictions 
of the SSF model and the measurements is reasonably good. The fact 
that turbulence levels roughly correspond to values estimated from a 
model which ignores effects of particle motion on turbulence prop- 
erties, indicates that effects of turbulence modulation (suggested 
by A1 Taweel and Landau [31]) or turbulence production by drops were 
small for the present flows. This is reasonable, since the compari- 
son is confined to dilute portions of the spray. 

Predicted and measured profiles of Reynolds stress are illus- 
treated in Fig. 22. The SSF predictions are adequate for both sprays. 
This is consistent with the reasonably good predictions of mean ve- 
locities and k obtained with this model for the test sprays. 

Effects of the presence of drops on turbulence properties are 
more evident when individual components of velocity fluctuations are 
examined. Measured radial profiles of u 1 , v 1 and w 1 are illustrated 
in Figs. 23 and 24 for the finely and coarsely atomized sprays. Pre- 
dictions were obtained assuming (u*2:v'2 ;w f ^) * (1:0. 5:0. 5) k, which 
is approximately observed in the fully-developed region of single- 
phase round jets [26,32]. Predictions constructed in this manner are 
in fair agreement with the measurements — particularly in the region 
far from the injector. A notable feature of the results, however, is 
that levels of anisotropy are rather high for positions near the 
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injector and generally exceed levels observed for comparable values 
of x/d in single-phase jets [26,32]. Since this region abuts the 
dense-spray portion of the flow, it seems likely that the presence 
of drops are responsible for the higher degree of anisotropy since 
effects of slip are preferentially transmitted into the streamwlse 
velocity component. High levels of anisotropy are also of concern 
regarding the prescription of eddy properties in the present SSF 
model, since this approach is based on the assumption of isotropic 
turbulence. This suggests that multistress models of particle-laden 
flovs might profitably be examined in order to gain more insight con- 
cerning effects of particles on turbulence properties. 


6. SUMMARY AND CONCLUSIONS 


New measurements of mean and fluctuating gas velocities, drop 
sizes and liquid fluxes were completed In the dilute portions of 
nonevaporating sprays in a still environment. These results were 
combined with existing measurements of the structure of particle- 
laden jets to evaluate models of the two-phase, jet-like flows. 

Three models were considered, as follows: (1) a locally homogeneous 
flow (LHF) model, where slip between the phases is neglected; (2) a 
deterministic separated flow (DSF) model, where interphase slip is 
considered but turbulent dispersion of particles is neglected; and 
(3) a stochastic separated flow (SSF) model, modified slightly from 
a proposal by Gosman and loannides [12], which allows for both inter- 
phase slip and the dispersion of particles by turbulent fluctuations. 
All three models employed a well-calibrated k-e model for predictions 
of continuous-phase properties [1,5-7,16]. Two versions of the SSF 
model were considered, one computing the transit time of a particle 
through an eddy using a linearized procedure similar Gosman and 
loannides [12], the other allowing for nonlinear effects in transit 
time evaluations. 

The LHF model was satisfactory for flows containing tracer-like 
particles, but was unsatisfactory for most practical flows. For 
particle-laden jets, the LHF model generally overestimated the rate 
of flow development and dispersion of particles, similar to past ex- 
perience with this model [5-7]. For the nonevaporating sprays, how- 
ever, the LHF model underestimated flow widths — unlike the results of 
earlier work. Such enhanced turbulent dispersion of particles, for 
certain ranges of turbulence and particle properties, is often ob- 
served in multiphase flows. The effect was more evident for the 
present sprays than for the particle-laden jets due to the smaller 
density of the liquid, which allowed the drops to respond more readily 
to turbulent fluctuations; and greater initial slip and rates of de- 
celeration in the sprays, due to smaller injector dimensions. Due to 
the possibility of preferential dispersion under some conditions, we 
conclude that the LHF model does not always provide an upper bound 
on the rate of development of sprays, as suggested in the past fl]. 
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The DSF model generally underestimated the rate of flow develop- 
ment and particle dispersion for the present data base. This approach 
appears to have limited utility for modeling practical particle-laden 
flows, Go8man and Ioannldes [12] note, however, that uncertainties in 
initial conditions for multiphase flows are potentially a greater 
source of error than neglecting particle dispersion* 

In contrast, the SSF models yielded reasonably good results over 
the present data base. In particular, the nonlinear version requires 
only a slightly more computational effort than the linearized version, 

while generally yielding improved predictions. The greatest dis- 
crepancies between predictions of the SSF models and the measurements 
in particle-laden jets were obseived at high particle loadings. Uncer- 
tainties in initial conditions of the experiments are a major factor 
In these discrepancies, e.g., potential Magnus effects discussed by 
Laats and Frishman [24,24] for their measurements and particle ve- 
locity fluctuations at the jet exit for the measurements of Levy and 
Lockwood [25], Effects of turbulence modulation, discussed by A1 Taweel 
and Landau [31] or turbulence generation by particles may also contribute 
to errors at high particle mass loadings. The SSF models also provided 
adequate treatment of enhanced particle dispersion in the spray with 
no modification of the model from its original calibration (where ef- 
fects of enhanced dispersion were not observed). While this is encour- 
aging, additional evaluation of the model is needed — particularly con- 
sidering improved specifications of initial conditions and measurements 
of drop sizes and velocities throughout the flow. Work along these 
lines is currently in progress in this laboratory. 

Present measurements were limited to the dilute portion of the 
flows, where void fractions were greater than 99.1%. In this region, 
major effects of particles/drops on turbulence properties were not 
observed. As the dense particulate flow region was approached, how- 
ever, turbulent velocity fluctuations exhibited increased anisotropy — 
suggesting a significant modification of turbulence properties by 
particles. The present SSF model, which employs assumptions of 
isotropy, still performed reasonably well in spite of this effect. 
However, further consideration of effects of anisotropy, perhaps 
using a multistress turbulence model, would be desirable. Further 
development of the model is also needed to consider direct effects of 
particles on the turbulence properties of the continuous phase. 

Based on experience to date, the SSF method appears to provide an 
attractive formulation for treating nonlinear interphase transport 
phenomena in sprays. Extension of this approach to treat processes 
of drop heat-up and evaporation in sprays is currently being considered 
in this laboratory. 
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APPENDIX A 


Data for the Nonevaporating Spray (Case 1) 


A,1 Gag Phase 


Table A.l Axial Variation of Centerline Velocity 
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u / u 
c o 
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50 
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400 
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Table A. 2 Radial Variation of Quantities 
x/d - 50 


r/x 

u/u 

c 

- 2 
U V7U 

c 

/Z^/Z 

c 

/v' 2 /u 

c 

/v’ 2 /u 

c 


0.0 

1.0 

0.0 

0.229 

0.104 

0.099 

.0365 

0.017 

0.983 

.0051 

0.231 

0.108 

0.101 

.0375 

0.034 

0.952 

.0091 

0.229 

0.110 

0.112 

.0386 

0.050 

0.886 

.0116 

0.238 

0.113 

0.118 

.0416 

0.067 

0.799 

.0134 

0.239 

0.124 

0.117 

.0431 

0.084 

0.722 

.0153 

0.238 

0.223 

0.112 

.0421 

0.101 

0.654 

.0132 

0.222 

0.1.16 

0.106 

.0371 

0.117 

0.557 

.0120 

0.221 

0.116 

0.099 

.0359 

0.134 

0.422 

.0104 

0.189 

0.109 

.189 

.0279 

0.151 

0.357 

.0082 

0.186 

0.108 

O'* 

.0264 

0.168 

0.277 

.0069 

0.178 

0.086 

0.069 

.0220 

0.184 

0.233 

.0060 

0.141 

0.080 

0.064 

.0151 

0.201 

0.180 

.0046 

0.133 

0.066 

0.057 

.0127 

0.218 

0.085 

.0025 

0.075 

0.060 

0.051 

.0059 
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Table A. 3 Radial Variation of Quantities 


x/d « 100 


r/x 

®/“ c 


/?^/u 

c 

/?*/u 

c 

/^/u 

c 

k/u 2 
c 

u'v’/u 1 
c 

0.0 

1.0 

0.0 

0.224 

0.147 

0.161 

.0489 

0.008 

0.988 

.0039 

0.224 

0.151 

0.163 

.0496 

0.017 

0.982 

.0053 

0.226 

0.151 

0.173 

.0521 

0.025 

0.961 

.0076 

0.236 

0.150 

0.165 

.0528 

0.034 

0.932 

.0100 

0.240 

0.151 

0.166 

.0541 

0.042 

0.885 

.0127 

0.240 

0.153 

0.169 

.0552 

0.050 

0.832 

.0157 

0.246 

0.156 

0.171 

.0571 

0.067 

0.717 

.0199 

0.247 

0.144 

0.174 

.0559 

0.084 

0.602 

.0202 

0.236 

0.146 

0.172 

.0533 

0.101 

0.515 

.0195 

0.228 

0.141 

0.165 

.0497 

0.117 

0.409 

.0175 

0.209 

0.137 

0.156 

.0433 

0.134 

0.310 

.0149 

0.186 

0.128 

0.138 

.0351 

0.151 

0.243 

.0119 

0.161 

0.116 

0.121 

.0270 

0.168 

0.168 

.0082 

0.142 

0.103 

0.101 

.0205 

0.184 

0.121 

.0068 

0.126 

0.090 

0.086 

.0157 

0.201 

0.059 

.0014 

0.061 

0.072 

0.053 

.0059 
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Table A. 4 Radial Variation of Quantities 
x/d « 250 


r/x 

u/u 

c 

— 7 — T/ — 2 

u’v'/u 

fz^tZ 

c 

/Z^/Z 

c 

/z^/z 

c 

k/u 2 
c 

0.0 

1.00 

0.0 

0.258 

0.207 

0.209 

.0764 

0.017 

0.981 

.0101 

0.257 

0.209 

0.211 

.0771 

0.034 

0.923 

.0146 

0.257 

0.210 

0.214 

.0781 

0.050 

0.799 

.0185 

0.257 

0.206 

0.216 

.0776 

0.067 

0.709 

.0202 

0.251 

0.207 

0.216 

.0762 

0.084 

0.597 

.0190 

0.239 

0.194 

0.212 

.0698 

0.101 

0.487 

.0189 

0.223 

0.185 

0.198 

.0616 

0.117 

0.389 

.0159 

0.200 

0.169 

0.184 

.0513 

0.134 

C.297 

.0129 

0.181 

0.152 

0.163 

.0411 

0.151 

0.191 

.0079 

0.147 

0.133 

0.142 

.0297 

0.168 

0.146 

.0055 

0.118 

0.109 

0.106 

.0185 

0.184 

0.104 

.0032 

0.091 

0.091 

0.078 

.0112 
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or root 



Table A. 5 Radial Variation of Quantitiea 


x/d ■ 600 


r/x 

u/u 

c 

u’v'/u ^ 
c 

/z^/z 

c 

/z^/z 

c 

« 

CM 

1^0° 

UA 

0.0 

1.00 

0.0 

0.272 

0.231 

.0906 

0.017 

0.997 

.0072 

0.277 

0.233 

.0928 

0.034 

0.957 

.0162 

0.281 

0.233 

.0936 

0.050 

0.869 

.0185 

0.274 

0.230 

.0905 

0.067 

0.778 

.0195 

0.269 

0.218 

.0841 

0.084 

0.666 

.0196 

0.259 

0.210 

.0777 

0.101 

0.546 

.0189 

0.237 

0.201 

.0684 

0.117 

0.446 

.0153 

0.219 

0.186 

.0589 

0.134 

0.313 

.0104 

0.187 

0.167 

.0454 

0.151 

0.228 

.0079 

0.163 

0.147 

.0349 

0.168 

0.166 

.0050 

0.116 

0.089 

.0148 

0.184 

0.087 

.0022 

0.070 

0.061 

.00612 
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A. 2 Liquid Phase 

Table A. 6 Axial Variation of Centerline Liquid Flux 


x/d 

G /G x 10 3 
c o 

40 

5.25 

50 

3.72 

70 

2.65 

100 

1.45 

150 

0.96 

250 

0.41 

400 

0.19 


600 


0.08 
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Table A. 7 Radial Variation of Liquid Flux 

r/x 

x/d - 50 

5/5 

c 

100 

G/G 

c 

250 

G/G 

c 

600 

G/G 

c 

0.0 

1.00 

1.00 

1.00 

1.00 

0.017 

— 

0.85 

0.97 

0.97 

0.034 

0.97 

0.80 

0.90 

0.90 

0.050 

— 

0.77 

0.75 

0.83 

0.067 

0.85 

0.75 

0.65 

0.67 

0.084 

— 

0.68 

0.44 

0.50 

0.101 

0.73 

0.56 

0.33 

0.32 

0.117 

— 

0.44 

0.20 

0.21 

0.134 

0.56 

0.38 

0.16 

0.12 

0.151 

— 

0.29 

0.12 

0.05 

0.168 

0.42 

0.21 

0.05 

— 

0.184 

— 

0.15 

0.01 

— 

0.201 

0.24 

0.08 

— 

— 

0.218 

— 

0.05 

— 

— 

0.235 

0.11 






0.268 


0.07 


Table A. 8 Radial Variation of Drop Size Distributions at x/d 
(Entries are Relative Percentage Number Frequency) 
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(Continued) 

Table A. 8 Radial Variation of Drop Size Distributions at x/d 
(Entries are Relative Percentage Number Frequency) 
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T^ble A. 9 Radial Variation of SMD at x/d • 50 


r/x 

SMD (ym) 

0.00 

40.6 

0.050 

41.8 

0.084 

46.2 

0.117 

50.7 

0.151 

49.3 

0.184 

9.4 

0.218 

43.6 


0.251 


47.6 
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APPENDIX B 


Data for the Nonevaporatiag Spray (Cage 2) 


B.l Gaa Phase 


Table B.l Axial Variation of Centerline Velocity 


x/d 

u /u 
c o 

40 

0.490 

50 

0.447 

70 

0.338 

100 

0.249 

150 

0.182 

250 

0.118 

400 

0.078 

500 

0.066 

600 

0.060 
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Table B.2 Radial Variation of Quantities 


x/d - 50 


r/x 

u/u 

c 

— r 2 
u f v f /u 

c 

/?*/g 

c 

/v’ 2 /u 

c 

/w' 2 /u 

c 

ir \ 

0.0 

1.000 

0.00 

0.228 

0.125 

0.126 

.0416 

0.017 

0.998 

.0034 

0.234 

0.127 

0.126 

.0434 

0.034 

0.956 

.0060 

0.236 

0,128 

0.127 

.0441 

0,050 

0.924 

.0103 

0.235 

0.132 

0.130 

.0448 

0.067 

0.829 

.0132 

0.222 

0.134 

0.132 

.0422 

0.084 

0.773 

.0152 

0.222 

0.136 

0.133 

.0429 

0.101 

0.693 

.0179 

0.226 

0.136 

0.133 

.0434 

0.117 

0.614 

.0196 

0.222 

0.134 

0.130 

.0420 

0.134 

0.550 

.0166 

0.220 

0.129 

0.128 

.0408 

0.151 

0.462 

.0140 

0.216 

0.127 

0.122 

.0388 

0.168 

0.391 

.0133 

0.198 

0.114 

0.105 

.0316 

0.184 

0.267 

.0080 

0.169 

0.104 

0.094 

.0241 

0.201 

0.148 

.0041 

0.110 

0.081 

0.078 

.0124 
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Table B.3 Radial Variation of Quantities 


x/d - 100 


r/x 

C 1 

C 1 
O 

u'v'/u * 
c 

/Pvg 

C 

Mg 

c 

c 

E/5 c 2 

0.0 

1.00 

0.00 

0.248 

0.157 

0.163 

.0564 

0.017 

0.998 

.0046 

0.247 

0.159 

0.159 

.0558 

0.034 

0.972 

.0096 

0.252 

0.161 

0.162 

.0576 

0.050 

0.899 

.0149 

0.258 

0.164 

0.166 

.0605 

0.067 

0.836 

.0179 

0.266 

0.166 

0.171 

.0637 

0.084 

0.745 

.0199 

0.269 

0.164 

0.172 

.0645 

0.101 

0.656 

.0192 

0.259 

0.162 

0.172 

.0613 

0.117 

0.548 

.0169 

0.252 

0.156 

0.166 

.0576 

0.134 

0.447 

.0142 

0.208 

0.149 

0.155 

.0447 

0.151 

0.333 

.0086 

0.181 

0.138 

0.141 

.0359 

0.168 

0.159 

.0054 

0.138 

0.117 

0.121 

.0237 

0.184 

0.113 

.0010 

0.089 

0.076 

0.088 

.0108 
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Table B.4 Radial Variation of Quantities 


x/d - 250 


r/x 

u/u. 

u'v'/u 2 
c 

Zz^/z 

c 

Zz^/Z 

c 

Zz^/z 

c 

1 N 

•5° 

9M 

0.0 

1.00 

0.0 

0.246 

0.192 

0.189 

.0663 

0.017 

0.983 

.0061 

0.248 

0.196 

0.197 

.0694 

0.034 

0.923 

.0141 

0.251 

0.193 

0.199 

.0700 

0.050 

0.863 

.0192 

0.253 

■ .195 

0.198 

.0706 

0.067 

0.739 

.0209 

0.251 

0.201 

0.203 

.0722 

0.084 

0.644 

.0203 

0.241 

0.203 

0.209 

.0716 

0.101 

0.517 

.0169 

0.234 

0.177 

0.203 

.0636 

0.117 

0.453 

.0149 

0.211 

0.169 

0.181 

.0530 

0.134 

0.336 

.0113 

0.176 

0.154 

0.167 

.0414 

0.151 

0.199 

.0067 

0.149 

0.146 

0.150 

.0331 

0.168 

0.121 

.0035 

0.099 

0.090 

0.101 

.0141 
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Table B.5 Radial Variation of Quantities 
x/d - 500 


r/x 

u/u 

c 

-r-r /- 2 
u'v'/u 

c 

/g?/j 

c 

c 

k/u 2 
c 

0.00 

1.000 

0.0 

0.255 

0.204 

. 0742 

0.017 

0.985 

.0091 

0.258 

0.202 

.0741 

0.034 

0.939 

.0148 

0.261 

0.207 

.0768 

0.050 

0.866 

.0194 

0.254 

0.205 

.0745 

0.067 

0.813 

.0206 

0.246 

0.187 

.0652 

0.084 

0.646 

.0180 

0.224 

0.180 

.0575 

0.101 

0.532 

.0160 

0.214 

0.159 

.0483 

0.117 

0.365 

.0118 

0.190 

0.140 

.0378 

0.134 

0.251 

.0089 

0.150 

0.119 

.0257 

0.151 

0.163 

.0063 

0.110 

0.088 

.0138 

0.168 

0.069 

.0024 

0.045 

0.046 

.0032 


Calculated assuming J w’ 2 - /v' 2 
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B.2 Liquid Phase 

Table B.6 Axial Variation of Centerline Liquid Flux 


x/d 

G /G x 10 3 
c o 

40 

3.24 

50 

2.13 

70 

1.34 

100 

0.99 

150 

0.50 

250 

0.25 

400 

0.12 

500 

0.076 
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Table B.7 Radial Variation of Liquid Flux 


r/x 

x/d - 50 
C/G c 

100 

G/G 

c 

250 

G/G 

c 

500 

G/C 

c 

0.0 

1.00 

1.00 

1.00 

1.00 

0.017 

— 

1.01 

0.99 

0.98 

0.034 

0.97 

0.99 

0.87 

0.92 

0.050 

— 

0.88 

0.80 

0.85 

0.067 

0.99 

0.81 

0.70 

0.80 

0.084 

— 

0.76 

0.52 

0.68 

0.101 

0.88 

0.63 

0.43 

0.55 

0.117 

— 

0.57 

0.32 

0.39 

0.134 

0.80 

0.45 

0.27 

0.27 

0.151 

— 

0.37 

0.20 

0.20 

0.168 

0.62 

0.31 

0.16 

0.12 

0.184 

— 

0.26 

0.14 

— 

0.201 

0.44 

0.20 

0.10 

— 

0.218 

— 

0.16 

0.07 

— 

0.235 

0.20 

— 

— 

— 

0.268 

0.13 





Table B.8 Radial Variation of Drop Size Distributions at x/d 
(Entries are Relative Percentage Number Frequency) 
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0.117 1.2 0.62 2.5 0.0 1.9 1.2 0.6 0.6 0.0 0.6 0.0 


Table B.8 Radial Variation of Drop Size Distributions at x/d 
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Table B.9 Radial Variation of SMD at x/d - 50 


r/x 

SMD (pm) 

0.00 

111.1 

0.050 

109.3 

0.084 

116.6 

0.117 

126.9 

0.151 

100.8 

0.184 

112.1 

0.218 

91.1 


